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ABSTRACT 

We present measurements of Fe, Mg, Si, Ca, and Ti abundances for 388 radial velocity member 
stars in the Sculptor dwarf spheroidal galaxy (dSph), a satellite of the Milky Way. This is the 
largest sample of individual a element (Mg, Si, Ca, Ti) abundance measurements in any single dSph. 
The measurements are made from Keck/DEIMOS medium-resolution spectra (6400-9000 A, i? ^ 
6500). Based on comparisons to published high-resolution {R > 20000) spectroscopic measurements, 
our measurements have uncertainties of crpe/H] = 0.14 and ^[a/Fe] — 0.13. The Sculptor [Fe/H] 
distribution has a mean ([Fe/H]) = —1.58 and is asymmetric with a long, metal-poor tail, indicative 
of a history of extended star formation. Sculptor has a larger fraction of stars with [Fe/H] < —2 than 
the Milky Way halo. We have discovered one star with [Fe/H] = — 3.80±0.28, which is the most metal- 
poor star known anywhere except the Milky Way halo, but high-resolution spectroscopy is needed to 
measure this star's detailed abundances. As has been previously reported based on high-resolution 
spectroscopy, [a/Fe] in Sculptor falls as [Fe/H] increases. The metal-rich stars ([Fe/H] ~ —1.5) have 
lower [a/Fe] than Galactic halo field stars of comparable metallicity. This indicates that star formation 
proceeded more gradually in Sculptor than in the Galactic halo. We also observe radial abundance 
gradients of —0.030 ± 0.003 dex per arcmin in [Fe/H] and -1-0.013 ± 0.003 dex per arcmin in [a/Fe] 
out to 11 arcmin (275 pc). Together, these measurements cast Sculptor and possibly other surviving 
dSphs as representative of the dwarf galaxies from which the metal-poor tail of the Galactic halo 
formed. 

Subject headings: galaxies: individual (Sculptor dwarf) — galaxies: dwarf — galaxies: abundances — 
Galaxy: evolution — Local Group 



1. INTRODUCTION 

The dwarf spheroidal galaxy (dSph) companions of the 
Milky Way (MW) are excellent laboratories for investi- 
gating the chemical evolution and star formation histo- 
ries of dwarf galaxies. These galaxies have undergone 

I at m ost a few star formation episodes (Holtzma n et a l. 

. I2006D and are dynamically simple fWalker et all 120071 ). 
The dSphs of the MW provide an opportunity to ex- 
amine closely the processes that establish the galaxy 
luminosity-metallicity relation fe.g.. lSalvadori fc Ferraral 
[2009). 

The MW dSphs are also considered to be strong candi- 
dates of a population of dwarf galaxies that were tidally 
stripped by the young Galaxy and eventually incorpo- 
rated into the Galactic halo. This scenario has be- 
come central to our p icture of how large ga laxies form 
(jSearle fc ZinnI 119781 : [Robertson et"al1 I20050 . Impor- 
tant tests of this scenario are to compare the details of 

^ Data herein were obtained at the W. M. Keck Observatory, 
which is operated as a scientific partnership among the California 
Institute of Technology, the University of California, and NASA. 
The Observatory was made possible by the generous financial sup- 
port of the W. M. Keck Foundation. 



the metallicity distribution function of the collection of 
dSphs to that of the Galactic halo stars and to compare 
abundance ratio patterns seen in dSphs to those mea- 
sured for the halo (e.g.. iVenn et al.ll20 04). 

To date, each of these areas has been hampered by 
the small sample of dSph stars for which high-quality 
measurements of [Fe/H] and ab undance ratios for other 
eleme nts have been available. iLanfranchi fc Matteuccil 
([2004*) compared their models of dSphs less massive than 
Sagittarius to six or fewer stars per galaxy. The usual 
approach for high-quality detailed abundance determina- 
tions is to use high-resolution spectroscopy (HRS, R > 
20000) of individual stars. Because of the large distances 
to even the nearest dSphs, these are time-consuming ob- 
servations even using the largest telescopes. 

Our approach is to derive abundances from medium- 
resolution spectroscopy (MRS, R ~ 6500) using the 
Deep Imaging M ulti-Object Spectrometer (DEIMOS, 
iFaber et al.ll2003D o n the Keck II telescope. As demon- 
strated bylKirbxei al. (2008a|l), accurate measurements 
can be made for Fe and some a elements (Mg, Si, Ca, and 
Ti) w ith these individual stellar spectra. Shetrone et al.l 
(|2009f ) demonstrated similarly precise results using the 
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Keck I LRIS spectrometer on a sample of individual stars 
in the Leo II dSph. In a typical dSph, the DEIMOS field 
of view allows between 80 and 150 red giant stars to 
be targeted per multi-object mask. Samples of several 
hundred giants can be observed in a given dSph. The 
Dwarf Abund ances and Radial Velocities team (DART, 
iTolstov et al.|[2b 04. hereafter TG4) has been collecting a 
combination of MRS and HRS in dSphs to exploit the 
advantages of both techniques. 

This paper is the first in a series that explores the 
multi-element abundances of stellar systems measured 
with MRS. The particular focus of this series is to char- 
acterize the distributions of [Fe/H] and [a/Fe] in MW 
dSphs. These measurements will provide insight into the 
role of dSphs in bu i lding the Galactic s tellar halo (i.e., 
ISearle fc ZinnHTflTl IWhite fc ReeslfTOTl) . 

Our first target is the S culptor dSph (a = I'^OO'", 
S = -33°43', My = -11.1. iMateol [19981 1. Sculptor has 
been a favored HRS and MRS target for the past ten 
years. Of all the dSphs, it appears most often in expla- 
nations of dSph chemical evolut ion and galaxy form ation 
(e.g., TG4, Sh etrone et al.l [2003':'Geis ler et al.l l2007^. T04 
discovered that Sculptor is actually "two galaxies" in one, 
with two stellar populatio ns that are kinematical ly and 
compositionally distinct. iBattaglia et al.l ()2006f ) later 
showed that Fornax also displays multiple stellar pop- 
ulations with different kinematics, spatial extents, and 
metallicities. But Sculptor is also uniq ue in that it is 
the only MW dS ph known to rotate (jBattaglia et al.l 
I2008al) . Recently. IWalker et all (|2009f l pub lished radial 
veloci t ies fo r 1365 Sculptor members, and iVenn fc Hilll 
(|2005l I2008D presented high-resolution abundance mea- 
surements of Mg, Ca, Ti, and Fe for 91 stars in Sculptor. 
They also measured Y, Ba, and Eu for some of those 
stars. 

This paper consists of six sections and an ap- 
pendix. Section [2] introduces the spectroscopic tar- 
get selection and observations, and Sec. [3] explains 
how the spectra are prepared for abundance mea- 
surements. Section |4] describes the technique to ex- 
tract abu ndances, which builds on the metho d de- 
scribed bv lKirbv. Guhathakurta. &z Sneded l|2008al here- 
after KGS08). In Sec. [5l we present the metalHcity dis- 
tribution and multi-element abundance trends of Sculp- 
tor. In Sec. [6l we summarize our findings in the con- 
text of dSph chemical evolution and the formation of the 
Galaxy. Finally, we devote the appendix to quantifying 
the uncertainties in our MRS measurements, including 
comparisons to independent HRS of the same stars. 

2. OBSERVATIONS 

2.1. Target Selection 

We s elected targets f rom t he Sculptor photometric cat- 
alog of iWestfall et al.l (|2006l ). The catalog includes pho- 
tometry in three filters: M and T2 in the Washing- 
ton system, and the intermediate-width DD051 filter 
(henceforth called D) centered at 5150 A. This band 
probes the flux from a spectral region susceptible to ab- 
sorpti on by the surfac e gravi ty-sei isitive Mg I an d MgH 
lines. iMaiewski et~an (|2000[ ) and IWestfall et all (|2006f ) 
outlined the procedure for distinguishing between distant 
red giant stars and foreground Galactic dwarf stars using 
these three filters. We followed the same procedure to se- 
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Fig. 1. — DEIMOS slitmask footprints laid over a map of sources 
from the photometric catalog. Targets selected for spectroscopy 
are shown in red. Targets observed in more than one mask are 
shown in green. Blue diamonds enclose stars with previous HRS 
abundance measurements. The left and bottom axis scales show 
the angular displacement in arcmin from the center of the galaxy 
(ao = l'^00™09=, <5o = -33°42'30", Mateo 1998), and the right 
and top axis scales s how the physical displac ement for an assumed 
distance of 85.9 kpc l lPietrzviiski et al.ll2008l) . 

lect a sample of red giant candidates from the Sculptor 
MT2D catalog. 

Nine stars, listed in Table [U have previously pub- 
lished HRS abundance measurements (jShetrone et al.l 
120031 iGeisler et all l2005f ). These stars were observed 
and provide the basis for demonstrating the accuracy 
of the MRS abundance measurements, described in the 
appendix. 

2.2. Slitmask Design 

We designed the DEIMOS slitmasks with the IRAF 
software module dsimulator.^ Each slitmask subtended 
approximately 16' x 4'. In order to adequately subtract 
night sky emission lines, we required a minimum slit 
length of 4". The minimum space between slits was 
0'.'35. When these constraints forced the selection of 
one among multiple possible red giant candidates, the 
brightest object was selected. The slits were designed 
to be at the approximate parallactic angle at the antici- 
pated time of observation (—25°). This choice minimized 
the small light losses due to differential atmospheric re- 
fraction. This configuration was especially important for 
Sculptor, which was visible from Keck Observatory only 
at a low elevation. The slitmasks' sky position angle (PA) 
was -35°. The 10° offset between the slit PA and the 
slitmask PA tilted the night sky emission lines relative to 
the CCD pixel grid to increase the subpixel wavelength 
sampling and improve sky subtraction. 

•^ [http: //www, ucolick. org/$~\slm$phlllips/deimos_ref /masks .html I 
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TABLE 1 

Targets with Previous High-Resolution Abundances 



iiciniG 


rGference 


RA 


Dec 


M 


To 


H482 


Shetrone et al. f2003") 


00''59™58=2 


-33°41'08" 


17.967 ±0.030 


16.324 ± 0.020 


H459 


Shetrone et al. f2003) 


Oihoomi2!5 


-33°43'01" 


18.465 ± 0.032 


16.924 ±0.031 


H479 


Shetrone et al. ('2003") 


0l'"00™12!7 


-33°41'15" 


17.562 ± 0.023 


15.860 ±0.030 


H400 


Shetrone et al. (2003) 


Oli^OO^lTfO 


-33°45'13" 


18.413 ± 0.030 


17.140 ±0.027 


H461 


Shetrone et al. f2003') 


0l'^00™18!2 


-33°42'12" 


17.806 ± 0.028 


16.166 ± 0.027 


1446 


Geisler et al. (2005) 


00'"59™46!4 


-33°41'23" 


17.618 ± 0.023 


15.695 ± 0.022 


195 


Geisler et al. (2005) 


00'i59™55?6 


-33°46'39" 


17.515 ± 0.022 


15.845 ± 0.018 


982 


&!islcr ct al. (2005) 


0l'^00™16!2 


-33°42'37" 


17.433 ± 0.025 


15.552 ± 0.028 


770 


Geisler et al. (2005 ~) 


0l'^00""23?8 


-33°42'17" 


17.623 ± 0.025 


15.857 ±0.026 



TABLE 2 
DEIMOS Observations 



Slitmask 


Targets 


UT Date 


Exposures 


Seeing 


sell 


86 


2008 Aug 3 


3 X 1200 s 


0'.'8 


scl2 


106 


2008 Aug 3 


2 X 900 s 


0'.'8 


scl3 


87 


2008 Aug 4 


1 X 462 s 


0"9 






2008 Aug 31 


1 X 1000 s 


0'.'8 






2008 Aug 31 


1 X 834 s 


0'.'8 


scl5 


95 


2008 Sep 1 


3 X 720 s 


0^8 


scl6 


91 


2008 Sep 1 


3 X 720 s 


1'.'2 



Note. — The scl4 slitmask was not observed. 



Figure [T] shows the coordinates of all the objects in 
the catalog regardless of their probability of membership 
in Sculptor. Five DEIMOS slitmask footprints enclose 
the spectroscopic targets: sell, scl2, scl3, scl5, and scl6 
(see Tab. [2]). The scl5 slitmask included 24 targets also 
included on other masks. These duplicate observations 
provide estimates of uncertainty in radial velocity and 
abundance measurements fSec. 13.3] and Sec. \A.l\ . The 
spectral coverage of each slit is not the same. The min- 
imum and maximum wavelengths of spectra of targets 
near the long, straight edge of the DEIMOS footprint can 
be up to 400 A lower than for targets near the irregularly 
shaped edge of the footprint (upper left and lower right 
of the slitmask footprints in Fig. [1] respectively). Fur- 
thermore, spectra of targets near either extreme of the 
long axis of the slitmask suffered from vignetting which 
reduced the spectral range. It is important to keep these 
differences of spectral range in mind when interpreting 
the differences of measurements derived from duplicate 
observations. 

Figure [2] shows the color-magnitude diagram (CMD) 
of the targets within the right ascension and declination 
ranges of the axes in Fig. [TJ The MT2D membership 
criteria caused the selected red giants to form a tight se- 
quence. This selection may have imposed a metallicity 
bias on the spectroscopic sample. Although only a tiny 
fraction of stars lay outside the main locus of the red 
giant branch, some may have been spectroscopically un- 
targeted members of Sculptor. For example, if Sculptor 
contained any old stars with [Fe/H] > —0.5, they would 
have been too red to be included in the spectroscopic 
sample. Any such metallicity bias should have excluded 
at most a few stars. 

2.3. Spectroscopic Configuration and Exposures 
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Fig. 2. — Color-magnitude diagram in the Washington and 
Cousins systems for the sources within the right ascension and 
declination ranges shown in Fig. [T] The symbols have the same 
meanings as in Fig. [T] The transformation from the Washington 
system (M and T2) to the Cousins syst em (Vc and Jp) is = T2 
and Vc~ Ic = 0.800(M - T2) - 0.006 l|Maiewski et al.|[200g ). 

Our observing strat egy was nearly identical to that of 
iSimon fc Gehal P OOTI and Kirbv^eTin (|2008aD . In sum- 
mary, we used with the 1200 lines mm~^ grating at a 
central wavelength of 7800 A. The slit widths were 0'.'7, 
yielding a spectral resolution of ^ 1.3 A FWHM (resolv- 
ing power R - 6500 at 8500 A). The OG550 fiher blocked 
diffraction orders higher than m = 1 . The spectral range 
was about 6400-9000 A with variation depending on the 
slit's location along the dispersion axis. Exposures of 
Kr, Ne, Ar, and Xe arc lamps provided wavelength cal- 
ibration, and exposures of a quartz lamp provided flat 
fielding. Table [2] lists the number of targets for each slit- 
mask, the dates of observations, the exposure times, and 
the approximate seeing. 

3. DATA REDUCTION 
3.1. Extraction of One- Dimensional Spectra 
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Fig. 3. — Examples of small regions of DEIMOS spectra of four 
different stars. The continuum in each spectrum has been normal- 
ized to unity. The Iq magnitude, measured effective temperature, 
and measured [Fe/H] is given for each star. The top two panels 
show two stars with very different [Fe /H] , and the bottom two pan- 
els show two stars with nearly the same temperature and [Fe/I^] 
but different SNR. The colors show the regions used to measure 
each of the Fe, Mg, Si, Ca, and Ti abundances (see Fig. |5j- 

We reduced the raw frames using version 1.1.4 of 
the DEIMOS data reduction pipeh ne developed by the 
DEEP Galaxy Redshift Survey.^ IGuhathakurta et all 
ffive the details of the data reduction. We also 
made use of th e optim izations to the code described by 
iSimon fc Gehal (|2007l . Sec. 2.2 of their article). These 
modifications provided better extraction of unresolved 
stellar sources. 

In summary, the pipeline traced the edges of slits in the 
flat field to determine the CCD location of each slit. The 
wavelength solution was given by a polynomial fit to the 
CCD pixel locations of arc lamp lines. Each exposure 
of stellar targets was rectified and then sky-subtracted 
based on a B-spline model of the night sky emission lines. 
Next, the exposures were combined with cosmic ray rejec- 
tion into one two-dimensional spectrum for each slit. Fi- 
nally, the one-dimensional stellar spectrum was extracted 
from a small spatial window encompassing the light of the 
star in the two-dimensional spectrum. The product of 
the pipeline was a wavelength-calibrated, sky-subtracted, 
cosmic ray-cleaned, one-dimensional spectrum for each 
target. 

Some of the spectra suffered from unrecoverable de- 
fects, such as a failure to find an acceptable polynomial 
fit to the wavelength solution. There were 53 such spec- 
tra. An additional 2 spectra had such poor signal-to- 
noise ratios (SNR) that abundance measurements were 
impossible, leaving 410 useful spectra, comprising 393 
unique targets and 17 duplicate measurements. 

Figure [3] shows four example spectra at a variety of 
Iq magnitudes, effective temperatures, and [Fe/H]. The 

^ |http : //astro .berkeley . edu/$~{\sim} $cooper/deep/spec2d/| 



two upper panels show stars in the top 10% of the SNR 
distribution. The two lower panels show stars from the 
middle and bottom 10% of the distribution. 

The one-dimensional DEIMOS spectra needed to be 
prepared for abundance measurements. The preparation 
included velocity measurement, removal of telluric ab- 
sorption, and continuum division. KGS08 (their Sec. 3) 
described these preparations in detail. We followed the 
same process with some notable exceptions, described 
below. 

3.2. Telluric Absorption Correction 

We removed the absorption introduced into the stellar 
spectra by the earth's atmosphere in the same manner 
as KGS08: division by a hot star template spectrum. 
However, the high airmass of the Sculptor observations 
caused much stronger absorption than KGS08 observed 
in globular cluster (GC) spectra. Even after scaling the 
hot star template spectrum by the airmass, large resid- 
uals in the Sculptor stellar spectra remained. Conse- 
quently, we masked spectral regions of heavy telluric ab- 
sorption before measuring abundances. These regions are 
6864-6932 A, 7162-7320 A, 7591-7703 A, 8128-8351 A, 
and 8938-10000 A (see Fig. [5]). 

3.3. Radial Velocities and Spectroscopic Membership 
Determination 

Our primary interest in this paper is chemical abun- 
dances, and we measured radial velocities only to deter- 
mine membership and to shift the spectra into the rest 
frame. 

Following KGS08, we measured stellar radial velocities 
by cross-correlation with a template spectrum. However, 
KGS08 cross-correlated the observed spectra against syn- 
thetic spectra whereas we cross-correlated the observed 
spectra against high SNR template spectra of stars ob- 
served with DEIMOS. Templates observed with the same 
instrument should provide more accurat e radial velocity 
measu rements than synthetic templates. ISimon fc Gehal 
()2007f ) provided their template spectra to us. For the 
rest of the analysis, the spectra are shifted to the rest 
frame. 

Although the MT2D selection eliminated almost all of 
the foreground MW contaminants from the spectroscopic 
sample, we checked the membership of each target by ra- 
dial velocity selection. Figure [3] shows the distribution of 
radial velocities in this spectroscopic data set along with 
the best-fit Gaussian. We consider the radial velocity 
limits of Sculptor membership to be 84.8 km < Vr < 
138.3 km s~^. We chose these limits because beyond 
them, the expected number of Sculptor members per 
2 km s^^ bin (the a pproxim ate maximum velocity res- 
olution of DEIMOS. [Simon fc Geha 2007 ) is fewer than 
0.5. This selection eliminated just 5 out of 393 unique 
targets. 

As a check on our procedure, we compared some de- 
rived quantities from the velocity distribution to previ- 
ous measurements. The mean velocity of our sample 
is (I'hciio) = 111.6 ± 0.5 km s~^ with a dispersion of 
(Ty = 8.0 ± 0.7 km s~^. The velocity dispersion is the 
per-measurement velocity error subtracted in quadrature 
from the la width of the velocity distribution. The per- 
measurement error is 3.9 km s~^, which is the standard 
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Fig. 4. — Distribution of measured radial velocities for targets 
in the Sculptor field along with the best-fit Gaussian. The top 
left label gives the mean and standard deviation of this Gaussian 
fit. The five stars outside of the dashed lines are not considered 
Sculptor members. The velocity range of this plot includes all stars 
for which a velocity measurement was possible. 

deviations of the differences in measure d velocities for 
the 17 duphcate spectra. In comparison, IWestfall et al.l 
((200l found (wheiio) = 110.4 ± 0.8 km s'^ (difference of 
+1.2(t) and cr„ = 8.8 ± 0.6 km s^^ (difference of -0.9cr). 
The comparison of the velocity dispersions depends on 
the assumed binary fraction (Queloz ct al. 1995) and — 
given the presence of multiple kinematically and spa- 
tially distinct populations in Sculptor (T04 ) — the region 
of spectrosc opic selection. Furthermore, I Walker et al.l 
([2007', '200!?) reported velocity dispersion gradients, and 
[Battaglia ct al. (2008a) reported mean velocity gradients 
along the major axis, indicating rotation. We choose not 
to address the kinematic complexity of this system in 
this paper. 

3.4. Continuum Determination 

In the abundance analysis described in Sec. Hj it is 
necessary to normalize each stellar spectrum by dividing 
by the slowly varying stellar continuum. KGS08 deter- 
mined the continuum by smoothing the regions of the 
stellar spectrum free from strong absorption lines. In- 
stead of smoothing, we fit a B-spline with a breakpoint 
spacing of 150 A to the same "continuum regions" de- 
fined by KGS08. Each pixel was weighted by its inverse 
variance in the fit. Furthermore, the fit was performed 
iteratively such that pixels that deviated from the fit by 
more than ba were removed from the next iteration of 
the fit. 

The spline fit results in a smoother continuum determi- 
nation than smoothing. Whereas the smoothed contin- 
uum value may be influenced heavily by one or a few pix- 
els within a relatively small smoothing kernel, the spline 
fit is a global fit. It is more likely to be representative of 
t he true stellar continuu m than a smoothed spectrum. 

iShetrone et al.l (|2009f) pointed out the importance of 
determining the continuum accurately when measur- 
ing weak lines in medium-resolution spectra. They re- 
fined their continuum determinations by iteratively fit- 



TABLE 3 

New Grid of ATLAS9 Model Atmospheres 



Parameter 


Minimum Value 


Maximum Value 


Step 


Tcff (K) 


3500 


5600 


100 


5600 


8000 


200 


logg (cm s~^) 


0.0 (Teff < 7000 K) 


5.0 


0.5 




0.5 (Teff > 7000 K) 


5.0 


0.5 


[A/H] 


-4.0 


0.0 


0.5 


[o/Fe] 


-0.8 


+1.2 


0.1 



ting a high-order spline to the quotient of the observed 
spectrum and the best-fitting synthetic spectrum. We 
adopted this procedure as well. As part of the iterative 
process described in Sec. 14.71 we fit a B-spline with a 
breakpoint spacing of 50 A to the observed spectrum di- 
vided by the best-fitting synthetic spectrum. We divided 
the observed spectrum by this spline before the next it- 
eration of abundance measurement. 

4. ABUNDANCE MEASUREMENTS 

The following section details some improvements on 
the abundance measurement techniques of KGS08. As- 
pects of the technique not mentioned here were un- 
changed from the technique of KGS08. In summary, 
each observed spectrum was compared to a large grid 
of synthetic spectra. The atmospheric abundances were 
adopted from the synthetic spectrum with the lowest x^. 

A major improvement was our measurement of four 
individual elemental abundances in addition to Fe: Mg, 
Si, Ca, and Ti. We chose these elements because they are 
important in characterizing the star formation history of 
a stellar population and because a significant number 
of lines represent each of them in the DEIMOS spectral 
range. 

4.1. Model Atmospheres 

Like KGS08, we built sy nthetic sp ectra based on AT- 
LAS9 model atmosp heres (iKuruc^ IT993f) with no con- 
vective overshooting (jCastelli et al.|[l997l ). KGS08 chose 
to allow the atmospheres to have [a/Fe] = -1-0.4 or 
[a/Fe] =0.0. This choice allowed them to use the large 
grid of ATLAS9 model atmosp heres computed with new 
opacity distribution functions (jCastelh fc Kuruc^l2004l ). 
However, we found that best-fitting model spectra com- 
puted by KGS08 tended to cluster around [a/Fe] = -1-0.2 
due to the discontinuity in caused by the abrupt 
switch between alpha-enhanced and solar-scaled models. 

To avoid this discontinuity, we recomputed ATLAS9 
model atmospheres on the grid summarized in Table [3l 
The new grid required recomputing new opacity distri- 
bution functions (ODFs), for which we used the DF- 
SYNTHE code (Castelll [20051 ). Unlike the grid of 
iCastelli fc Kuruczr(l2004D. we adopt ed the solar compo- 
sition of 'Anders fc Gr evessi (I1989D. except for Fe, for 
which we followed Sncd en et al.l()1992 . see the note in Ta- 
ble [4]). One opacity distribution function was computed 
for each of the 189 combinations of [A/H] and [a/Fe] 
specified in Table [H The abundances of all the elements 
except H and He were augmented by [A/H]. Addition- 
ally, the abundances of O, Ne, Mg, Si, Ar, Ca, and Ti 
were augmented by [a/Fe]. These ODFs were used to 
compute one ATLAS9 model atmosphere for each grid 
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point in Table [3] and for two values of microturbulent 
velocity, for a total of 139104 model atmospheres. 

4.2. Microturbulent Velocity 

In order to reduce the number of parameters required 
to determine a stellar abundance, KGS08 assumed that 
the microturbulent velocity (^) of the stellar atmosphere 
was tied to the surface gravity (logg). They chose to fit 
a line to the spe ctroscopica ll y mea sured ^ and log 5 of 
the giant stars in iFulbrightl 's (|200ClD sample: 

C (km s"^) = 2.70-0.51 logg (1) 

We also adopted a relation between ^ and \ogg, but we 
re-determined this relatio n from t h e GC r ed giant sample 
of KGS08 combined with iKirbvT s ()2009D compilation of 
high-resolution spectroscopic measurements from the lit- 
erature (Frcbcl ct al. 2009; Geisler et al. 2005; Johnson 
[2?)0l (Lai et al.ll2007l : IShetrone et al.ll2n0ll . 1200.1 and ref- 
erences from KGS08). The best-fit line between the spec- 
troscopically measured ^ and log 5 is 

^ (km s"^) = (2.13 ± 0.05) - (0.23 ± 0.03) logg . (2) 

corresponding roughly to a 0.0-0.5 km s^^ decrease in 
^, depending on logg. In the generation of the grid of 
synthetic stellar spectra described in Sec. 14.41 ^ was not 
a free parameter, but was fixed to logg via Eq. [2l 

In general, a decrease in ^ increases the measurement 
of [Fe/H]. Therefore, this change tended to increase the 
derived values of [Fe/H]. A typical change in [Fe/H] was 
< -1-0.05 dex. This change would be more severe in an 
HRS analysis based on equivalent widths (EWs). In our 
minimization, the abundance measurement was most 
sensitive to lines with large d(EW) / d[Fe/ll]. Such lines 
are the weak, unsaturated transitions whose strength 
does not depend on ^. The DEIMOS spectra contain 
enough of these weak lines that did not play a large 
role in the abundance determination. 

4.3. Line List 

We compared the Fe I oscillator strengths (logg/) in 
the KGS08 line list to va lues measured in the labora- 
tory ()Fuhr fc 'wresell2006( ). Most of the KGS08 oscilla- 
tor strengths were stronger than the laboratory measure- 
ments. The average offset was 0.13 dex. Because KGS08 
calibrated their line list to the solar spectrum, we in- 
terpreted this offset as a systematic error in the solar 
model atmosphere, solar spectral synthesis, and/or solar 
composition. Accepting the laboratory-measured values 
as more accurate than the so lar calibration, we replaced 
Fe I oscillator strengths with iFuhr fc Wiesd where avail- 
able, and we subtracted 0.13 dex from logg/ for all other 
Fe I transitions in the KGS08 line list. All other data re- 
mained unchanged. 

Decreasing the oscillator strengths requires a larger 
[Fe/H] to match the observed spectrum. The amount of 
change in [Fe/H] depends on the atmospheric parameters 
as well as the saturation of the measured Fe lines. From 
comparison of results with the old and new line lists, we 
estimate a typical change in [Fe/H] to be ~ -1-0.1 dex. 

4.4. Generation of Synthetic Spectra 



TABLE 4 

Adopted Solar, Composition 



Element 


12 + loge 


Element 


12 + log e 


Mg 


7.58 


Ti 


4.99 


Ca 


6.36 


Fe 


7.52 


Si 


7.55 







Note. — This co mposition is adopted from 
lAnders Hz Grev csse ( 1989|) , except for Fe. For 
justificatio n of the adopted F e solar abun- 
dance, see ISneden et all l(T99l ). The abun- 
dance of an element X is defined as its number 
density relative to hydrogen: 12 + log ex = 
12 -I- log(nx) - log(nH). 

The spectra were synthesized as described in KGS08. 
Specifically, the current version of the local thermody- 
namic e quilibrium (L TE) spectrum synthesis software 
MOOG (|Snedenlll973( ) generated one spectrum for each 
point on the grid. The spectral grid was more finely 
spaced in [Fe/H] than the model atmosphere grid. The 
spacing is 0.1 dex for each of [Fe/H] and [a/Fe], yielding 
a total of 316848 synthetic spectra. 

The solar composition used in the generation of the 
synthetic spectra was identical to the solar composition 
used in the computation of the model atmospheres. Ta- 
ble [4] lists the adopted solar abundances for the five el- 
ements for which we measure abundances in Sculptor 
stars. 

4.5. Effective Temperatures and Surface Gravities 

Different spectroscopic studies of chemical abundances 
rely on different sources of information for determin- 
ing the effective temperature (Toff) and surface grav- 
ity (logg) of the stellar atmosphere . KGS08 con- 
sulte d Yonsei-Yale model isochrones (jDemarque et al.l 
l2004f ) to determine the temperature and gravity that 
correspond to a dereddened color and an extinction- 
corrected absolute magnitu de. They al so considered 
Victoria- Regina (VandcnBe rg et ahl |2006l) and Padova 
(jGirardi et al.l l2002) model isochrones, as well as an em- 
pirica l color-temperature relation (Ramirez fc Melended 
[2005h . 

The Fe lines accessible in DEIMOS spectra span a large 
range of excitation potential. Together, these different 
lines provide a constraint on Toff. KGS08 (their Sec. 5.1) 
showed that — without any photometric information — 
the synthesis analysis of medium-resolution spectra of 
GC stars yielded values of Toff very close to values previ- 
ously measured from HRS. Therefore, we chose to mea- 
sure Toff from photometry and spectroscopy simultane- 
ously. 

To begin, we converted extinction-corrected 
(jSchlegel et al.lll998h Washington M an d T2 magnitudes 
to C ousins Vc and Ic magnitudes (jMaiewski et al.l 
[2000h . With these magnitudes, we computed Tefi from 
the Yonsei-Yale, Victoria- Regina, and Padova mode l 
isochrones, as well as the iRamfrez fc MelendezI ()2005[ ) 
empirical color-based Toff. For each measurement, we es- 
timated the effect of photometric error by measuring the 
standard deviation of Tgff determined from 1000 Monte 
Carlo realizations of Vc and Ic- In each realization, Vc 
and Ic were chosen from a normal distribution with a 
mean of the measured, extinction-corrected magnitude 
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and a standard deviation of the photometric error. We 
call this error STeff^i, where i represents each of the four 
photometric methods of determining Toff. In order to 
arrive at a single photometric Tcff, we averaged the four 
Tcff,i together with appropriate error weighting. We also 
estimated the random and systematic components of 
error. In summary, 



I ^i-e = 



^cff — ■ 



E.ST J, 



(3) 

(4) 



(^totalToff = ^ ((^randTcff)^ + ((^sysTcff)^ 



(6) 



For the stars in this data set, the median random, 
systematic, and total errors on Toff were 98 K, 58 K, 
and 117 K respectively. The somewhat large errors on 
the photometric temperatures indicated that the spec- 
tra may help constrain Toff. Therefore, Eq. [3] does not 
show the final temperature used in the abundance deter- 
mination. Section [4771 describes the iterative process for 
determining Toff and elemental abundances from spec- 
troscopy. 

We followed a similar procedure for determining logg 
photometrically, except that we used only the three 
model isochrones and not any empirical calibration. 
The error on the true distance modulus (19.67 ± 0.12, 
iPietrzvnski et al1l2008f ) was included in the Monte Carlo 
determination of the error on logg. The median ran- 
dom, systematic, and total errors on log 5 were 0.06, 0.01, 
and 0.06. These errors are very small, and the medium- 
resolution, red spectra have little power to help constrain 
log g because there are so few ionized lines visible. There- 
fore, we assumed the photometric value of log g for the 
abundance analysis. 

4.6. Wavelength Masks 

The procedure described in the next section consisted 
of separately measuring the abundances of five elements: 
Mg, Si, Ca, Ti, and Fe. The procedure relied on finding 
the synthetic spectrum that best matched an observed 
spectrum. In order to make this matching most sensitive 
to a particular element, we masked all spectral regions 
that were not significantly affected by abundance changes 
of that element. 

To make the wavelength masks, we began with a base 
spectrum that represented the solar composition in which 
the abundances of all the metals were scaled down by 
1.5 dex ([A/H] = —1.5). The temperature and gravity 
of the synthetic star were Toff — 4000 K and logg — 1.0. 
Then, we created two pairs of spectra for each of the 
five elements. In one spectrum, the abundance of the 
element was enhanced by 0.3 dex, and in the other, de- 
pleted by 0.3 dex. Spectral regions where the flux differ- 
ence between these two spectra exceeds 0.5% were used 
in the abundance determination of that element. This 
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Fig. 5. — A coaddition of all Sculptor stars with the continuum 
normalized to unity. The high SNR provided by the coaddition 
makes stellar absorption lines readily apparent. The colored re- 
gions show the wavelength masks used in the determination of the 
abundance of each element. Regions susceptible to telluric absorp- 
tion are labeled with blue text. Because large residuals from the 
telluric absorption correction remain, we eliminate these regions 
from the abundance analysis. Some stellar features excluded from 
the abundances measurement are also labeled. 

small threshold assured that weak lines, which experi- 
ence large fractional changes in EW as [Fe/H] changes, 
were included in the analysis. We repeated this proce- 
dure for spectra with Toff = 5000 K, 6000 K, 7000 K, and 
8000 K. Additional spectral regions that passed the 0.5% 
flux difference criterion were also included in the abun- 
dance determination of that clement. All other wave- 
lengths were masked. 

The result was one wavelength mask for each of Mg, 
Si, Ca, Ti, and Fe, shown in Fig. [51 We also created 
one "a" mask as the intersection of the Mg, Si, Ca, and 
Ti masks. The a element regions do not overlap with 
each other, but the a element regions do overlap with 
the Fe regions. The most severe case is the Ca mask, 
where ~ 35% of the pixels are shared with the Fe mask. 
However, the overlap did not introduce interdependence 
in the abundance measurements. The a element abun- 
dances were held fixed while [Fe/H] was measured, and 
the Fe abundance was held fixed while [a/Fe] was mea- 
sured. The measurements of [Fe/H] and [a/Fe] were per- 
formed iteratively (see the next subsection). We tested 
the independence of the measurements by removing all 
overlapping pixels from consideration. Abundance mea- 
surements changed on average by only 0.01 dex. 

4.7. Measuring Atmospheric Parameters and Elemental 
Abundances 

A Levenberg-Mar quardt algorithm (the IDL routine 
MPFIT, written bv iMarkward^ 120091 ) found the best- 
fitting synthetic spectrum in ten iterative steps. In each 
step, the was computed between an observed spec- 
trum and a synthetic spectrum degraded to match the 
resolution of the observed spectrum. First, we interpo- 
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lated the synthetic spectrum onto the same wavelength 
array as the observed spectrum. Then, we smoothed 
the synthetic spectrum through a Gaussian filter whose 
width was the observed spectrum's measured resolution 
as a function of wavelength. 

1. Tes and [Fe/H], first pass: An observed spectrum 
was compared to a synthetic spectrum with Tcff 
and log 5 determined as described in Sec. 14.51 and 
[Fe/H] determined from Yonsei-Yale isochrones. 
For this iteration, [a/Fe] was fixed at 0.0 (solar), 
and only spectral regions most susceptible to Fe ab- 
sorption (Sec. l4.6() were considered. The two quan- 
tities Tcff and [Fe/H] were varied, and the algo- 
rithm found the best-fitting synthetic spectrum by 
minimizing x^- We sampled the parameter space 
between grid points by linearly interpolating the 
synthetic spectra at the neighboring grid points. 
Teff was also loosely constrained by photometry. As 
the spectrum caused T^s to stray from the photo- 
metric values, increased, and it increased more 
sharply for smaller photometric errors (as calcu- 
lated in Eq. [6]). Therefore, both photometry and 
spectroscopy determined Tcff. Photometry alone 
determined logg. 



a/Fe], first pass: For this iteration, Toff, logg, and 
Fe/H] were fixed. Only [a/Fe] was allowed to vary. 

In the model stellar atmosphere, the abundances of 

the a elements with respect to Fe varied together. 

Only the spectral regions susceptible to absorption 

by Mg, Si, Ca, or Ti were considered. 



3. Continuum refinement: The continuum-divided, 
observed spectrum was divided by the synthetic 
spectrum with the parameters determined in steps 
1 and 2. The result approximated a flat noise spec- 
trum. To better determine the continuum, we fit 
a B-spline with a breakpoint spacing of 50 A to 
the residual spectrum. We divided the observed 
spectrum by the spline fit. 

4. [Fe/H], second pass: We repeated step 1 with the 
revised spectrum, but Tcs was held fixed at the 
previously determined value. 

5. [Mg/Fe]: We repeated step 2. However, only Mg 
spectral lines were considered in the abundance 
measurement. 

6. [Si/Fe]: We repeated step 5 for Si instead of Mg. 

7. [Ca/Fe]: We repeated step 5 for Ca instead of Mg. 

8. [Ti/Fe]: We repeated step 5 for Ti instead of Mg. 

9. [a/Fe], second pass: We repeated step 2 for all 
of the a elements instead of just Mg. This step 
was simply a different way to average the a el- 
ement abundances than combining the individual 
measurements of [Mg/Fe], [Si/Fe], [Ca/Fe], and 
[Ti/Fe]. 

10. [Fe/H], third pass: The value of [a/Fe] affected the 
measurement of [Fe/H] because [a/Fe] can affect 



TABLE 5 

Systematic Abundance Errors 



Element Ratio 


"^sys 


Element Ratio 






Fe/H] 


0.136 


[Ca/Fe] 


0.087 




Mg/Fe] 


0.108 


[Ti/Fe] 


0.101 




Si/Fe] 


0.179 







the structure of the stellar atmosphere. Specif- 
ically, the greater availability of electron donors 
with an increased [a/Fe] ratio allows for a higher 
density of H~ ions. The subsequent increase in con- 
tinuous opacity decreases the strength of Fe and 
other non-a element lines. With [a/Fe] fixed at 
the value determined in step 9, we re-measured 
[Fe/H]. Typically, [Fe/H] changed from the value 
determined in step 1 by much less than 0.1 dex. 

4.8. Correction to [Fe/H] 

In comparing our MRS measurements of [Fe/H] to HRS 
measurements of the same stars (see the appendix), we 
noticed that our measurements of metal-poor stars were 
consistently ' -^ 0.15 dex lower. T he same pattern is also 
visible in the lKirbv et"aLl (|2008aD GC measurements (see 
their Figs. 6, 7, 10, and 11). 

We have thoroughly examined possible sources of this 
difference of scale. The changes to the microturbulent 
velocity relation (Sec. 14. 2p and the line list (Sec. 14. 3p 
were intended to yield a more accurate and standardized 
estimation of [Fe/H], but the offset still remained. Re- 
stricting the analysis to narrow spectral regions did not 
reveal any systematic trend of [Fe/H] with wavelength. 

A possible ex planation for this offset is overionization 
iTh evenin &: Id iart 1999). Ultraviolet radiation in stellar 
atmospheres can ionize Fe more than would be expected 
in LTE. Therefore, the abundance of Fe I would seem to 
be lower than the abundance of Fc II in an LTE analysis. 
Fe II does not suffer from this effect. However, the effect 
is smaller at higher [Fe/H], and we do not observe a trend 
with metallicity for the offset of our values relative to 
HRS studies. 

In order to standardize our measurements with previ- 
ous HRS studies, we added 0.15 dex to all of our mea- 
surements of [Fe/H]. This offset and the microturbu- 
lent velocity-surface gravity relation are the only ways in 
which previous HRS studies inform our measurements. 
Furthermore, this offset is not intended to change the 
standardization of our abundances. All of the abundance 
in this article, including those from other studies, are 
given relative to the solar abundances quoted in Table SI 

4.9. Error Estimation 

We repeated the error estimation procedure described 
by KGS08 (their Sec. 6) by repeating their abundance 
analysis on GC stars with the above modifications. We 
no longer found a convincing trend of (5[Fc/H] with 
[Fe/H]. Instead, we estimate the total error on [Fe/H] by 
adding a systematic error in quadrature with the SNR- 
dependent uncertainty of the synthetic spectral fit. The 
magnitude of (5sys[Fe/H] = 0.136 was the value required 
to force HRS and MRS [Fe/H] estimates of the same GC 
stars to agree at the la level. We also estimated sys- 
tematic errors for each of [Mg/Fe], [Si/Fe], [Ca/Fe], and 
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Fig. 6. — The metallicity distribution in Sculptor. The red curve 
is the maximum likehhood fit to a galactic chemical evolution 
model with pre-enrichment (Eq.[Tll, and the green curve is the max- 
imum likelihood fit to a model of star formation in the presence of 
infalling, zero-metallicity gas (Eq. [TTJ. The long, metal-poor tail 
is typical for systems with non-instantaneous star formation. 

[Ti/Fe] in tlie same manner as for [Fe/H]. These are listed 
in Table [3 

5. RESULTS 

In this section, we discuss the interpretation of the 
abundance measurements in Sculptor, all of which are 
presented in Table [6] on the last page of this manuscript. 

5.1. Metallicity Distribution 

The metallicity distribution function (MDF) of a dwarf 
galaxy can reveal much about its star formation his- 
tory. In chemical evolution models of dwarf galax- 
ies (e . g., [L anfranchi & Matteucci 2004; Marcolini et ahj 
120061 12008 '') ■ the duration of star formation affects the 
shape of the MDF. The MDF also has implications for 
the for mation of the MW. If the MW halo was b uilt from 
dSphs (|Searle fc ZinnlflQTl IWhite fc ReeslllQTl . then it 
is important to find dSph counterpar ts to halo field stars 
at all metallicities, as pointed out bv lHelmi et al.l (|2006l 
hereafter HOG). 

Figure [5] shows the MDF of Sculptor. The shape of the 
MDF is highly asymmetric, with a long, m etal-poor tail 
(as predicted bv lSalvadori fc Ferrarall2009l ) . The inverse- 
variance weighted mean is ([Fe/H]) = —1.58 with a stan- 
dard deviation of 0.41. The median is —1.58 with a me- 
dian absolute deviation of 0.33 and an interquartile range 
of 0.67. 

The MDF boasts an exceptionally metal-poor star, 
S1020549. The metallicity is [Fe/H] = -3.80 ± 0.28. 
Figure [7] shows how weak the Fe absorption lines are in 
this star. Frebel, Kirby, fc Simon (in preparation) have 
confirmed this extremely low metallicity with a high- 
resolution spectrum. 

Sculptor is now the most luminous dSph in which 
an extremely metal -poor (EMP , [Fe/H ] < —3) star 
has been detected. [Kirbv et al.l ([2008b) discovered 15 
EMP stars acro s s eigh t ultra- faint dwarf galaxies, and 
iCohen fc Huand (|2009( ) discovered one EMP star in the 
Draco dSph.] Stars more metal-poor than S1020549 are 
known to exist only in the field of the Milky Way field. 
This discovery hints that dSph galaxies like Sculptor may 
have contributed to the formation of the metal-poor com- 
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Fig. 7. — Regions of the DEIMOS spectrum of the extremely 
metal-poor star S1020549, which has [Fe/H] = -3.80 ± 0.28. The 
spectrum appears particularly noisy because the y-axis range is 
small. Some Fe absorption lines are barely detectable, but all to- 
gether, they contain enough signal to make a quantitative mea- 
surement of [Fe/H]. The shading corresponds to the same spectral 
region shown in Fig. [5] Frebel, Kirby, & Simon (in preparation) 
will present a high-resolution spectrum of this star, which confirms 
the extremely low metallicity. 

ponent of the halo. We discuss Sculptor's link to the halo 
further in Sec. 15.1.31 

The MT2D photometric selection of spectroscopic tar- 
gets may have introduced a tiny [Fe/H] bias. Figure [2] 
shows that the RGB is sharply defined in Sculptor. Be- 
cause the number density of stars redward and blueward 
of the RGB is much lower than the number density on the 
RGB, the number of very young or very metal-poor stars 
(blueward) or very metal-rich stars (redward) missed by 
photometric pre-selection must be negligible. Further- 
more, the hard color cut (as opposed to one that depends 
onM -D color) was 0.6 < (M - T2)o < 2.2. The CMD 
gives no reason to suspect Sculptor RGB members out- 
side of these limits, but it is possible that some extremely 
blue Sculptor members have been excluded. 

5.1.1. Possible Explanation of the Discrepancy with 
Previous Results 

Our measured MDF and our detection of EMP stars in 
Sculptor are at odds with the findings of H06. Whereas 
our MDF peaks at [Fe/H] ^ —1.3, theirs peaks at 
[Fe/H] ^ —1.8. Furthermore, our observed MDF is much 
more asymmetric than that of HOG, which may even be 
slightly asymmetric in the opposite sense (a longer metal- 
rich tail). The greater symmetry would indicate a less 
extended star f ormation histor y or early infall of a large 
a mount of gas (i Prantzos 2003!). 

iBattaglia et al . (2008b, hereafter B08b) observed a 
subset of the HOG stars at high resolution. The MDFs 
from the two studies have noticeably different shapes. 
Figure [5] shows that the HRS MDF peaks at [Fe/H] - 
— 1.3, which is also the peak that we observe. The mean 
and standard deviation of their MDF are —1.56 and 0.38. 
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Fig. 8. — Sculptor's metallicity distribution as observed in this 
study (MRS, black) and by B08b (HRS, red), which is a subset 
of the MDF observed by H06 (Ca triplet, green). The CaT-based 
MDF is more metal-poor probably because the sample of H06 is 
more spatially extended than the other two samples. 

However, the MDF of H06 peaks at [Fe/H] 1.8, and 

the mean and standard deviation are —1.82 and 0.35. 
The overlapping stars between the samples of B08b and 
H06 agree very well. 

The most likely explanation for the different MDFs is 
the different spatial sampling of the three st udies. Sculp- 
tor h as a steep radia l metallicit y gradient (Tolstov et al.l 
[200l iWestfall eFall [2006 : .Walker et all JOOa, also see 
Sec. 15. 3p . The stars in the center of Sculptor are more 
metal-rich than stars far from the center. HOG s ampled 
stars out to the tidal radius (rt = 76.5 arcmin, iMateol 
IT99I . but we and B08b sampled stars only out to about 
11 arcmin. As a result, the mean metallicity of the H06 
CaT sample is lower than our MRS sample and the B08b 
HRS sample. In the next subsection, we address the 
chemical evolution of Sculptor based on its MDF. Our 
conclusions are based only on stars within the central 
11 arcmin. 

5.1.2. Quantifying Chemical Evolution in Sculptor 

In chemical evolution models, ext ended sta r forma tion 
produces a long, metal-poor tail. iPrantzoi (|2008n de- 
scribed the shape of the differential metallicity distribu- 
tion derived from a "simple model" of galactic chemical 
evolution. Expressed in terms of [Fe/H] instead of metal 
fraction Z, the predicted distribution is 

d[Fe/H] V ^ H P / 

(7) 

where p is the effective yield in units of the solar metal 
fraction {Zq) and [Fe/H]^ is the initial gas metallicity. 
An initial metallic ity is nee ded to resolve the Galactic 
G dwarf problem (Ivan den~B crgh 1962; Schmidt 196^). 
A is a normalization that depends on p, [Fe/H]^, the 
final metallicity [Fe/H]j, and the number of stars in the 
sample N: 

, ^ {NlnlO)/p 

exp(-M^^)-exp(-Mi!^) 

The red curve in Figure [6] is the two-parameter, maxi- 



mum likelihood fit to Eq. \7\ The likelihood Li that star 
i is drawn from the probability distribution defined by 
Eq. [7| is the integral of the product of the error distri- 
bution for the star and the probability distribution. The 
total likelihood L = Yli Li- The most likely p and [Fe/H]g 
are the values that maximize L. For display, the curve 
has been convolved with an error distribution, which is 
a composite of N unit Gaussians. N is the total number 
of stars in the observed distribution, and the width of 
the i"^ Gaussian is the estimated total [Fe/H] error on 
the i**^ star. This convolution approximates the effect 
of measurement error on the model curve under the as- 
sumption that the error on [Fe/H] does not depend on 
[Fe/H]. This assumption seems to be valid because our 
estimates of (5[Fe/H] do not show a trend with [Fe/H]. 

The most likely yield — largely determined by the 
[Fe/H] at the peak of th e MDF— i s y = O.OSIZq. 
[From the MDF of H06, iPrantzosI (j2008l ) calculated 
p = 0.016Zq.] We also measure [Fe/H]^ = -2.92. 
H06 also measured [Fe/H]„ = -2.90 ± 0.21 for Sculp- 
tor, even though they included stars out to the tidal ra- 
dius, which are more metal-poor on average than the 
centrally concentrated stars in our sample. (Instead of 
finding the maximum likelihood model, they performed 
a least-squares fit to the cumulative metallicity distri- 
bution without accounting for experimental uncertainty. 
In general, observational errors exaggerate the extrema 
of the metallicity distribution, and the least-squares fit 
converges on a lower [Fc/HJq than the maximum likeli- 
hood fit.) One explanation that they proposed for this 
non-zero initial metallicity was pre-enrichment of the in- 
terstellar gas that formed the first stars. Pre-enrichment 
could result from a relatively late epoch of formation 
for Sculptor, after the supernova (SN) ejecta from other 
galaxies enriched the intergalactic medium from which 
Sculptor formed. However, our observation of a star at 
[Fe/H] = —3.80 is inconsistent with pre-enrichment at 
t he level of [ Fe/H1„ = -2.9. 

iPrantzoi (|2008[ ) instead interpreted the apparent 
dear th of EMP sta rs as an indication of early gas in- 
fall (|Prantzosll2003f l. wherein star formation begins from 
a small amount of gas while the majority of gas that will 
eventually form dSph stars is still falling in. In order to 
test this alternative to pre-enrichment, we h ave also fit an 
Infall Model, the "Best Accretion Model" of lLvnden-Belil 
IT975L also see IPagel' 1997). It is one of the models which 
accounts for a time-decaying gas infall that has an ana- 
lytic solution. The model assumes that the gas mass g 
in units of the initial mass is related quadratically to the 
stellar mass s in units of the initial mass: 




where M is a parameter greater than 1. When A/ = 1, 
Eq. [9] reduces to 5 = 1 — s, which describes the Closed 
Box Model. Otherwise, M monotonically increases with 
the amount of gas infall a nd with the departu re fro m the 
Simpl e Model. Following iLvnden-Belil (|1975f) and iPagell 
(|1997f ). we assume that the initial and infalling gas metal- 
licity is zero. The differential metallicity distribution is 
described by two equations. 
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Equation [TU] is transcendental, and it must be solved for 
s numerically. Equation [TT] decouples the peak of the 
MDF from the yield p. As M increases, the MDF peak 
decreases independently of p. 

The green line in Fig. [6] shows the most likely Infall 
Model convolved with the error distribution as described 
above. The Infall Model has M = 1.76, which is only a 
small departure from the Simple Model. 

Neither the Simple Model nor the Infall Model fits 
the data particularly well. Both models fail to repro- 
duce the sharp peak at [Fe/H] ^ —1.3 and the steep 
metal-rich tail. However, the Infall Model does repro- 
duce the metal-poor tail about as well as the Simple 
Model. Therefore, the Infall Model is a reasonable al- 
ternative to pre-enrichment, and it allows the existence 
of the star at [Fe/H] = —3.80. In reality, a precise ex- 
planation of the MDF will likely incorporate the radial 
metallicity gradients and multiple, superposed popula- 
tions. It is tempting to conclude from Fig. [6] that Sculp- 
tor displays two metallicity populations. We have not 
attempted a two-component fit, but that would seem to 
be a re asonable appr oa ch for future work, especially in 
light of lTolstov et all 's (|2004f) report of two distinct stel- 
lar populations in Sculptor. 

Searches for the lowest metallicity stars in the MW 
halo have revealed so me exquisitely me tal-poor stars 
(e.g., [Fe/H] = -5.96. IWebel et al.ll2fl?i8h . Such exotic 
stars have not yet been discovered in any dSph. How- 
ever, if Sculptor was not pre-enriched, a large enough 
sample of [Fe/H] measurements in Sculptor — and possi- 
bly other dSphs — may reveal stars as metal-poor as the 
lowest metallicity stars in the MW halo. 

5.1.3. Comparison to the Milky Way Halo MDF 

ISearle fc Zin^ (fl97ll and IWhite fc Reei ([l978h 
posited that the MW halo formed from the accretion 
and dissolution of dwarf galaxies. The dSphs that 
exist today may be the survivors fr om the can nibalistic 
construction of the Galactic halo. Hel mi et al.l (j200^ 
suggested that at least some of the halo field stars could 
not have come from counterparts to the surviving dSphs 
because the halo field contained extre mely metal-poor 
stars w hereas the dSphs do not. However. iSchoerck et al.l 
(|2008') showed that the Hamburg/ESO Survey's halo 
MDF, after correction for selection bias, actually looks 
remarkably like the MDFs of the dSphs Fornax, Ursa 
Minor, and Draco. Furthermore, iKirbv et al.l (|2008br ) 
presented MRS evidence for a large fr action of EMP 
stars in the ultra-faint dSph sample of iSimon fc Gehal 
(|2b07i) . suggesting that today's surviving dSphs contain 
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Fig. 9. — The metal-poor tails o f the MDFs in Sculp tor {black) 
and Galactic halo field stars {red, lijchoerck et al.ll2008) shown as 
cumulative distributions, all normalized to the number of stars 
with [Fe/H] < —2. The green line shows the MDF measured by 
the DART team (Helmi et al. 2006) with a calibration based on 
the Ca triplet. The calibration may overpredict very low metallic- 
ities. The synthesis-based metallicities {black, this work) are valid 
at lower [Fe/H] than the Ca triplet [Fe/H]. Regardless, the halo 
has a steeper metal-poor tail than Sculptor in both representa- 
tions. Galaxies such as Sculptor were probably not the dominant 
contributors to the halo. 

stars that span the full range of metallicities displayed 
by the Galactic field halo population. 

We revisit the halo comparison with the present 
MDF for Sculptor. Figure [9] shows the metal-poor 
tail ([Fe/H] < -2) of the MRS synthesis-based Sculp- 
tor JVIDF presente d her e, the CaT-based Sculptor 
MDF (Helmi et al.. [20061) . and the MW halo MDF 
(jSchoerck et all 12008^ As observed in the com par- 
isons to other dSphs presented by ISchoerck et al.l the 
halo seems to have a steeper metal-poor tail than the 
CaT-based Sculptor MDF, despite the evidence that 
CaT-based metallicities overpredict [Fe/H] at [F e/H] < 
-2.2 fe.g- iKoch et all 120081 ? INorris etaLl 12001) . The 
synthesis-based MDF does not rely on empirical calibra- 
tions, and the techniq ue has been shown t o work at least 
down to [Fe/H] = -3 (|Kirbv et al.ll2008bD . 

This MDF shows that the halo has a much steeper 
metal-poor tail than Sculptor. This result is consistent 
with a merging scenario wherein several dwarf galax- 
ies significantly larger than Scul ptor contributed most 
of the stars to th e halo field (e.g.. iRobertson et al.ll2005l : 
iFont et al.l [20061 ) ■ In these models, the more luminous 
galaxies have higher mean metallicities. Galaxies with 
a Sculptor-like stellar mass are minority contributors to 
the halo field star population. Less lu minous galaxies 
are even more metal- poor (Kirby e t al l [2008b) . There- 
fore, Sculptor conforms to the luminosity-metallicity re- 
lation for dSphs, and the difference between Sculptor's 
MDF and the MW halo MDF does not pose a problem 
for hierarchical assembly. 

5.2. Alpha Element Abundances 

The discrepancy between halo and dSph abundances 
extends beyond the MDF. In the first HRS study 
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Fig. 10. — Multi-element abundances in Sculptor (black). The 
point sizes reflect the quadrature sum of the errors on [Fe/H] and 
[X/Fe], where larger points have smaller errors. The bottom panel 
shows the average of the four elements shown in the other panels. 
For comparison, the red error bars show the means and standard 
deviations from the seven GCs of KGS08. Because the Sculptor 
and GC abundances were measured in the same way, the compar- 
ison demonstrates that [X/Fe] declines with increasing [Fe/H] in 
Sculptor, but not the GCs. 



of stars in a dSph, iShetrone. Bolte fc Stetson! ()1998f ) 
found that the [Ca/Fe] ratio of metal-poor stars in 
Draco appeared solar, in contras t to th e enha nced 
halo field stars . IShetrone. Cote & Sargent! ()2001[ ) and 
IShetrone et"aI1 (|2003n confirmed the same result in Sex- 
tans, Ursa Minor, Sculptor, Fornax, Carina, and Leo I, 
and they included other a elements in addition to Ca. 

Here, we present the largest sample of [a/Fe] measure- 
ments in any dSph. Figure [TOl shows [Mg/Fe], [Ca/Fe], 
and [Ti/Fe] versus [Fe/H] for Sculptor. The figure also 
shows the mean and standard deviations of all of the in- 
dividual stellar abundance measurements for each of the 
seven GCs in the sample of KGS08. All of the modifi- 
cations to the KGS08 technique described in Sees. [3] and 
[4] apply to the GC measurements in Fig. [TOl Although 
our discussion in the appendix demonstrates that our 
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Fig. 11. — As in Fig. 1101 the black points show medium-resolution 
multi-element abundances in Sculptor. Th e points with erro r bars 
show published high-resolution data (Sh etrone et al.l 120031 . blue, 
and Geisler et al. 200^, green). The green line is the inverse 
variance-weighted average of at least 20 stars within a window of 
A[Fe/Hl = 0.25. The red line shows the chemical evolution model 
of Lanfranchi & Matteucci (2004, updated 2009). The onset of 
Type la SNe causes the decline in [X/Fe] with [Fe/H]. Mg declines 
steadily because it is produced exclusively in Type II SNe, but Si, 
Ca, and Ti are produced in both Type la and II SNe. 

measurements are accurate on an absolute scale by com- 
paring to several different HRS studies in Sculptor, it is 
also instructive to compare abundances measured with 
the same technique in two types of stellar systems. All 
four element ratios slope downward with [Fe/H] in Sculp- 
tor but remain flat in the GCs. Additionally, the larger 
spread of [Mg/Fc] than other clement ratios in the GCs is 
not due to larger measurement uncertainties but to the 
known intrinsic s pread of Mg abunda nce in some GCs 
(see the review bv lGratton et al.|[20iDl . [Si/Fe], [Ca/Fe], 
and [Ti/Fe] are more slightly sloped than [Mg/Fe] in 
Sculptor because both Type la and Type II SNe pro- 
duce Si, Ca, and Ti, but Typ e II SNe are almost solely 
responsible for producing Mg (jWooslev fc Weaver IPOS'). 
Finally, to maximize the SNR of the element ratio mea- 
surements, we average the four ratios together into one 
number called [a/Fe]. The [a/Fe] ratio is flat across the 
GCs, but it decreases with increasing [Fe/H] in Sculptor. 
Quantitative models of chemical evolution in dwarf 
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Fig. 12. — As in Fig. 1101 the black points show medium-resolution 
multi-element abundances in Sculptor. The colored po ints show 
different components of the Milky Way llVenn et al.l200^ : the thin 
disk (red), the thick disk (green), and the field halo (cyan). The 
dashed lines are moving averages of the MW data in 0.75 dex bins 
of [Fe/H]. In Sculptor, [a/Fe] falls at lower [Fe/H] than in the halo, 
indicating that the halo field stars were less polluted by Type la 
SNe and therefore formed more rapidly than Sculptor stars. 

galaxies are consistent with these trends. At a certain 
time corresponding to a certain [Fe/H] in the evolution 
of the dSph, Type la begin to pollute the interstellar 
medium with gas at subsolar [a/Fe]. JVIore metal-rich 
stars that form from this gas will have lower [g/Fe] 
than t he more metal-poor stars. iLanfranchi fc Matteuccil 
(|2004f ) have developed a sophisticated model that in- 
cludes SN feedback and winds. They predicted the 
abundance distributions of six dSphs, including Sculp- 
tor. Figure [TT] shows our measurements with their pre- 
dictions (updated with new SN yields, G. Lanfranchi 
2009, private communication). As predicted, the range 
of [JVIg/Fe] is larger than the range of [Ca/Fe] or [Si/Fe] 
because Mg is produced exclusively in Type II SNe 
whereas Si and Ca are produced in both Type la and 
II SNe (jWooslev fc Weaved IT995I) . We do not observe 
strong evidence for a predicted sharp steepening in slope 
of both elements at [Fe/H] ~ —1.8, but observational 
errors and intrinsic scatter may obscure this "knee." 
Also, the observed [Fe/H] at which [IVIg/Fe] begins to 
drop is higher than the model predicts, indicating a 



less intense wind than used in the model. Note that 
the element ratios [X/Fe] become negative (subsolar) 
at high enough [Fe/H] , as p redicted by the models. 
ILanfranchi fc MatteuccI ()2004f l do not predict [Ti/Fe] 
because it behaves more like an Fe-peak element than 
an a element. 

In addition to tr ends of [a/Fe] with [Fe/H], 

llMarcohni et~aLl (|2006l . 120081) predicted the distribution 
functions of [Fe/H] and [a/Fe] of a Draco-like dSph. The 
range of [Fe/H] they predicted is nearly identical to the 
range we observe in Sculptor, and the shapes of both 
distributions are similar. The outcome of the models de- 
pends on the mass of the d Sph. Sculpt or is ten times 
more luminous than Draco (|Mateolfl998f ) and therefore 
may have a larger total mass. [However, iStrigari et al.l 
((2008) find that all dSphs have the same dynamical mass 
within 300 pc of their centers. It is unclear whether the 
total masses of the original, unstripped dark matter ha- 
los are the same.] In principle, these chemical evolution 
models could be used to measure the time elapsed since 
different epochs of star formation and their durations. 
We defer such an analysis until the advent of a model 
based on a Sculptor- like luminosity or mass. 

In Fig. 1121 we compare individual stellar abundances 
in Sculptor to MW h alo and disk field stars (compilation 
bv lVenn et all l2004f) . As has been seen m many pre- 
vious studies of individual stellar abundances in dSphs, 
[a/Fe] falls at a significantly lower [Fe/H] in Sculptor 
than in the JVIW halo. The drop is particularly appar- 
ent in [JVIg/Fe], which is the element ratio most sensitive 
to the ratio of the contributions of Type II to Type la 
SNe. The other element ratios also drop sooner in Sculp- 
tor than in the halo, but appear lower than in the halo 
at all metallicities. Along with the JVIDF comparison in 
Sec.l5.1.3[ this result is consistent with the suggestion by 
[Robertson et al.l ()2005f ) that galaxies significantly more 
massive than Sculptor built the inner MW halo. Their 
greater masses allowed them to retain more gas and expe- 
rience more vigorous star formation. By the time Type la 
SNe diluted [a/Fe] in the massive halo progenitors, the 
metallicity of the star-forming gas was already as high 
as [Fe/H] = -0.5. In Sculptor, the interstellar [Fe/H] 
reached only —1.5 before the onset of Type la SNe pol- 
lution. 

5.3. Radial Abundance Distributions 

Because dSphs interact with the MW, t hey can lose 
gas th rough tidal or ram pressure stripping (jLin fc Fabeil 
Il983f ). The gas preferentially leaves from the dSph's out- 
skirts, where the gravitational potential is shallow. If 
the dSph experiences subsequent star formation, it must 
occur in the inner regions where gas remains. Sculp- 
tor's MDF suggests a history of extended star formation. 
Sculptor might then be expected to exhibit a radial abun- 
dance gradient in the sense that the inner parts of the 
dSph are more metal-rich than the outer parts. 

The detection of a radial metallicity gradient in 
Sculptor has been elusive. In a photometric study, 
[Hurley-Keller (2000) found no evidence for an age or 
metallicity gradient. Based on HRS observat i ons o f 
five stars (the sain e sample as IShetrone et al.l I2003D . 
iTolstov et al.l ()2003f ) found no correlation between [Fe/H] 
and spatial position. Finally, in a sample of 308 stars 
with CaT-based metallicities, T04 detected a significant 
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Fig. 13. — Spatial abundance distributions in Sculptor. Point 
sizes are larger for stars with smaller measurement uncertainties. 
The red points reflect the mean values in 1 arcmin bins, along with 
the errors on the means. The red lines are the least-squares linear 
fits. We detect a gradient of —0.030 it 0.003 dex per arcmin in 
[Fe/H] and +0.013 ± 0.003 dex per arcmin in [a/Fc]. 

segregation in Sculptor: a centrally concentrated, rela- 
tively metal-ricli component an d an ex tended, relatively 
metal-poor component. West fa ll et al. (f2006i) arrived at 
the same conclusion, and Walk er et al I ()2009f ) confirmed 
the existence of a [Fe/H] gradient in a sample of 1365 
Sculptor members. 

In order to detect a gradient, those studies targeted 
Sculptor stars at distances of more than 20 arcmin. The 
maximum elliptical radius of this study is 11 arcmin. 
Therefore, this study is not ideally designed to detect 
radial gradients. Figure fT3l shows the radial distribution 
of [Fe/H] and [a/Fe] in Sculptor. The x-axis is the length 
of the semi-major axis of an ellipse de fined b y Sculptor's 
position angle and ellipticity ( Mateol [199 8*1. Although 
this study is limited in the spatial extent of targets, we do 
detect a gradient of —0.030 ±0.003 dex per arcmin. This 
estimate is very clo se to the gradient observed by T04. 
IWalker et all (l2009h measure a shallower gradient, but 
they present their results against circular radius instead 
of elliptic al radius . 

Marcoli ni et al.l ([20 OS^ predict radial gradients in both 
[Fe/H] and [a/Fe] in dSphs. In particular, they expect 
shallower [Fe/H] gradients for longer durations of star 
formation. The gradient we observe is stronger than any 
of their models. They also expect very few stars with 
low [a/Fe] at large radius. Given that [a/Fe] decreases 
with [Fe/H] and [Fe/H] decreases with distance, it seems 
reasonable to expect that [a/Fe] increases with radius. In 
fact, we detect an [a/Fe] gradient of -1-0.013 ± 0.003 dex 
per arcmin. 



6. CONCLUSIONS 

Sculptor is one of the best-studied dwarf spheroidal 
satellites of the Milky Way. In the past ten years, at least 
five spectroscopic campaigns at both low and high res- 
olution have targeted this galaxy. More than any other 
dSph, Sculptor has aided in the understanding of the 
chemical evolution of dSphs and the construction of the 
Milky Way stellar halo. 

We have sought to increase the sample of multi-element 
abundances in Sculptor through MRS. The advantages 
over HRS include higher throughput per resolution ele- 
ment, the ability to target fainter stars, and multiplex- 
ing. The large sample sizes will enable detailed compar- 
isons to chemical evolution models of [a/Fe] and [Fc/H] 
in dSphs. The disadvantages include larger uncertain- 
ties, particularly for elements with few absorption lines 
in the red, and the inability to measure many elements 
accessible to HRS. MRS is not likely to soon provide in- 
sight into the evolution of neutron-capture elements in 
dSphs. 

In order to make the most accurate measurements pos- 
sible, we hav e made a number of improvemen t s to th e 
technique of iKirbv. Guhathakurta. fc SnedenI ()2008aD . 
We have consulted independent HRS of the same stars 
to confirm the accuracy of our measurements of [Fe/H], 
[Mg/Fe], [Ca/Fe], and [Ti/Fe]. In the case of [Fc/H] 
and the average [a/Fe] our MRS measurements are only 
slightly more uncertain than HRS measurements. 

Some of the products of this study include 

1. An unbiased metallicity distribution for 
Sculptor. Because the synthesis-based abun- 
dances do not rely on any empirical calibration, 
their applicability is unrestricted with regard to 
[Fe/H] range. The MDF is asymmetric with a long, 
metal-poor tail, as predicted by chemical evolution 
models of dSphs. Furthermore, fits to simple chem- 
ical evolution models shows that Sculptor's MDF 
is consistent with a model that requires no pre- 
enrichment. 

2. The largest sample of [a/Fe] and [Fe/H] 
measurements in any single dSph: 388 stars. 

We have confirmed the tr end for [a/Fe] to dec rease 
with [Fe/H], as shown bv lGeisler et all (l2007f) with 
just nine st ars from the studies of iShetrone et all 
(2003) and iGeisler eTall p005). Chemical evolu- 
tion models may be constructed from these mea- 
surements to quantify the star formation history of 
Sculptor. 

3. The detection of radial [Fe/H] and [a/Fe] 
gradients. Our sample probes a smaller range 
than previous studies; nonetheless, we find a 
—0.030 ± 0.003 dex per arcmin gradient in [Fe/H] 
and a -1-0.013 ± 0.003 dex per arcmin gradient in 
[a/Fe]. 

4. The discovery of a Sculptor member star 
with [Fe/H]= -3.80 ± 0.28. This discovery sug- 
gests that since-disrupted galaxies similar to Sculp- 
tor may have played a role in the formation of the 
Milky Way metal-poor halo. High-resolution spec- 
troscopy of individual stars will confirm or refute 
this indication. 
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Fig. 14. — Cumulative distribution of differences between the re- 
peat measurements of [Fe/H] for 17 stars divided by the estimated 
error of the difference. The curve is the integral of a unit Gaussian. 
The curve matches the distribution well, indicating that the errors 
are estimated properly. 

Much more can be done with this technique in 
other galaxies. The stellar population of a dSph 
depends heavily on its stellar mass. For instance, 
iLanfranchi &: Matteuccil (j2004f) and [Robertson et al.l 
(|2005( ) predict that more massive satellites have an [a/Fe] 
"knee" at higher [Fe/H]. In the next papers in this series, 
we intend to explore the multi-element abundance distri- 
butions of other dSphs and compare them to each other. 
We will observe how the shapes of the MDFs and the 
[a/Fe]-[Fe/H] diagrams change with dSph luminosity or 
stellar mass. These observations should aid our under- 
standing of star formation, chemical evolution, and the 
construction of the Galaxy. 
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Fig. 15. — Same as Fig. 1141 for [a/Fe], which is the average of 
[Mg/Fe], [Si/Fe], [Ca/Fe], and [Ti/Fe]. 
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APPENDIX 

ACCURACY OF THE ABUNDANCE MEASUREMENTS 

In order to quantify the accuracy of the MRS measurements, we examine the spectra of stars observed more than 
once and stars with previous HRS measurements. 



Duplicate Observations 

The repeat observations of 17 stars provide insight on the effect of random error on the measurements of [Fc/H] 
and [a/Fe]. Figures 114! and [T51 summarize the comparisons of measurements of different spectra of the same stars. 
They show the cumulative distribution of the absolute difference between the measured [Fe/H] and [a/Fe] for each 
pair of spectra divided by the expected error of the difference (see Sec. 14.91) . The solid curve is the integral of a unit 
Gaussian, which represents the expected cumulative distribution if the estimated errors accurately represent the true 
measurement errors. In calculating the expected error of the difference, we apply the systematic error to only one of 
the two stars. Even though the same technique is used to measure abundances in both stars, some systematic error is 
appropriate because the wavelength range within a pair of spectra differs by 300-400 A. The different Fe lines in these 
ranges span a different range of excitation potentials, and the Levenberg-Marquardt algorithm converges on different 
solutions. 



Comparison to High- Resolution Measurements 

The most reliable test of the MRS atmospheric parameter and abundance estimates is to compare with completely 
independent obser vations and analyses of the same stars. Table |S] lists the previous HRS measurements of nine 
Sculptor members (iShetro ne et al.ll2003l : iGeisler et al.ll2005l ) as well as the DEIMOS measurements of the same stars. 
Unfortunately, these two HRS studies share no stars in common and therefore cannot be compared with each other. 
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TABLE 6 

Abundances of Stars with Previous High- Resolution Spectroscopy 



star 




log 3 


5 


[Fo/H] 


[Mg/Fo] 


[Si/Fe] 


[Ca/Fe] 


[Ti/Fe] 




(K) 


(cm s~^) 


(km s 


1) (dex) 


(dex) 


(dex) 


(dex) 


(dex) 




Previous High-Resolution Measurements 


H482 


4400 


1.10 


1.70 


-1.24 ±0.07 


+0.09 ± 0.14 


-0.07 ±0.17 


+0.06 + 0.06 


-0.17 ±0.09 


H459 


4500 


1.00 


1.65 


-1.66 ±0.07 


+0.36 ± 0.14 


+0.22 + 0.17 


+0.24 + 0.05 


-0.05 ±0.09 


H479 


4325 


0.70 


1.70 


-1.77 ±0.07 


+0.26 ± 0.17 


+0.00 + 0.23 


+0.17 + 0.05 


-0.05 ±0.09 


H400 


4650 


0.90 


1.70 


-1.98 ±0.08 


+0.37 + 0.13 




+0.38 + 0.09 


-0.07 ±0.15 


H461 


4500 


1.20 


1.70 


-1.56 ±0.07 


+0.18 ± 0.12 


+0.14 + 0.18 


+0.22 + 0.06 


+0.00 + 0.09 


1446 


3900 


0.00 


2.30 


-1.22 ±0.13 


-0.17 ±0.15 


-0.33 + 0.18 


-0.12 + 0.17 


-0.14 + 0.15 


195 


4250 


0.20 


1.80 


-2.12 ±0.15 


+0.71 ± 0.15 


+0.03 + 0.18 


+0.14 + 0.13 


+0.12 + 0.15 


982 


4025 


0.50 


2.20 


-0.99 ±0.16 


-0.20 ±0.11 




-0.14 + 0.17 


-0.12 + 0.19 


770 


4075 


0.00 


1.90 


-1.74 ±0.13 


+0.18 ± 0.15 


-0.31 + 0.18 


+0.19 ± 0.13 


+0.19 + 0.15 


Medium-Resolution Measurements 


H482 


4347 


0.83 


1.95 


-1.24 ±0.14 


-0.03 ± 0.17 


-0.16 ±0.20 


-0.19 ± 0.13 


-0.12 + 0.11 


H459 


4390 


1.12 


1.88 


-1.88 ±0.14 




-0.41 ±0.48 


+0.34 + 0.27 


+0.08 + 0.18 


H479 


4271 


0.63 


1.99 


-1.79 ±0.14 


+0.27 + 0.42 


-0.22 ±0.22 


+0.24 + 0.18 


-0.13 + 0.11 


H400 


4692 


1.36 


1.82 


-1.97 ±0.15 




+0.58 + 0.23 


+0.90 + 0.65 


+0.10 + 0.19 


H461 


4313 


0.78 


1.96 


-1.81 ±0.15 


+0.36 + 0.61 


+0.05 + 0.32 


+0.39 ± 0.32 


+0.25 + 0.16 


1446 


3838 


0.49 


2.03 


-1.22 ±0.14 


-0.03 ± 0.20 


-0.07 + 0.21 


-0.05 ± 0.24 


-0.40 + 0.11 


195 


4308 


0.65 


1.99 


-2.05 ±0.14 


+0.41 ± 0.17 


+0.06 + 0.19 


+0.09 ± 0.10 


+0.02 + 0.11 


982 


4147 


0.52 


2.02 


-0.84 ±0.14 






-0.30 ±0.23 


-0.47 + 0.11 


770 


4247 


0.59 


2.00 


-1.62 ±0.14 


+0.13 ± 0.29 


-0.23 + 0.21 


+0.02 ± 0.11 


-0.04 + 0.11 



Of Toff, logg, and ^, any spectroscopic ab undance measu r ement is most sensitive to Teff. In general, underestimating 
Teff leads to an underestimate of [Fe/H]. iShetrone et al.l ()2003l . hereafter S03) deter mine Teg spec t rosco pically by 
minimizing the slope of the derived abundance for each line versus excitation potential. Gcisle r et al.l (|2005l . hereafter 
G05) determine Toff photometrically with empirical color-temperature relations. Figure [16] shows Teff from those 
studies and this one for each of the nine stars in common. The MRS temperatures do not follow the temperatures of 
either HRS study better than the other. 

Both 803 and G05 measure log;; spectroscopically from demanding ionization equilibrium: [Fc/H] measured from 
Fe I lines must match that measured from Fe II lines. However, our red spectra have very few measurable Fe II lines. 
Alternatively, logy may be determined from a star's absolute magnitude and Toff via the Stefan-Boltzmann law. Even 
though gravity depends on the inverse square of Teff and the inverse square root of luminosity, luminosity imposes 
a stronger constraint on logg because of its larger range on the RGB than Teff. Even accounting for the error in 
the distance modulus to Sculptor, the typical error on photometric logg is ~ 0.1 dex. Therefore, we determine logg 
from photometry alone. Figure [17] shows the comparison between logg used by S03 and G05 and this study. The 
agreement is not particularly good, with discrepancies up to 0.6 dex. However, the photometric values of logg are 
more accurate than can be determined from the medium-resolution red spectra, which show very few lines of ionized 
species. Furthermore, as discussed below, errors in logy influence the abundance measurements much less than errors 
in Teff. 

Both S03 and G05 measure microturbulent velocity (^) by forcing all Fe lines to give the same abundance regardless 
of their reduced width. We have fixed ^ to log g with an empirical relation (Eq. ^ . Figure [18] compares the HRS 
microturbulent velocities (^) to our adopted values. The largest discrepancy is 0.3 km s~^. 

Figure [19] shows the comparison between HRS and MRS [Fe/H] measurements for the same stars. The agreement is 
very good {a = 0.14 dex). Just two stars out of nine do not fall within la of the one-to-one line. 

The MRS [Fe/H] for star 770 is larger than the HRS [Fe/H]. The MRS T^e is also significantly larger than the 
HRS Teff for this star. Similarly, the MRS T^s for star H461 is lower than the HRS Teff, forcing the MRS [Fe/H] 
lower than the HRS [Fe/H]. In fact, even the smaller deviations from the [Fe/H] one-to-one line can be attributed to 
deviations from the Teff one-to-one line. No such correlation can be attributed to deviations in logy or ^. The close 
correspondence between Figs. 1 16l and 1 1 91 demonstrates that Teff is the dominant atmospheric parameter in determining 
metallicity. 

B08b pubhshed a catalog of VLT/FLAMES [Fe/H] measurements based on both the EW of the infrared Ca II triplet 
(CaT) and HRS (Hill et al., in preparation). The two resolution modes of FLAMES {R - 6500 and R ~ 20000) 
allowed them to complete both MRS and HRS analyses with the same instrument. Their high-resolution spectroscopic 
sample and ours overlap by 47 stars, which are shown in Fig. [20] The agreement {a = 0.14 dex) is as good as the 
previous comparison to HRS studies. 

The BOSb HRS measurements rely on atmospheric parameters determined from both five-band photometry and 
spectroscopy. We also measure Teff spectrophotometrically. Our methods may be similar, although we do not use 
infrared photometry. There appears to be a small systematic trend such that our MRS measurements are lower than 
the BOSb HRS measurements of [Fe/H] at both low and high [Fe/H]. The average discrepancy at the extrema of the 
residuals is 0.2 dex. We withhold a detailed investigation of these residuals until publication of the details of the HRS 
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Fig. 16. — Comparison between effective temperature (T^g) used 
in previous HRS abundance analyses and photometric T^f{ used for 
this work's MRS abundance analysis. Symbol shape indicates the 
reference for the HRS abundances. Star names from Tab. [l] are 
printed to the upper left of each point. 



Fig. 17. — Same as Fig. [16] except for surface gravity (logg). The 
error bars represent photometric error, and they are given by replac- 
ing Tcff with log g in Eq. \E\ 



study. 

B08b share seven stars in common with SOS and G05. To emphasize the accuracy of our MRS analysis, we note that 
the scatter of the differences between the two sets of HRS studies {a = 0.16 dex) is in fact larger than the scatter in 
the comparison between the MRS [Fe/H] and the same seven stars of S03 and G05 {a = 0.14 dex). This small sample 
does not indicate that the MRS measurements are more accurate than any HRS measurements, but it does suggest 
that the accuracy is competitive. 

Figure!^ shows the comparison between the MRS and HRS (SOS and G05) values of [Mg/Fe], [Si/Fe], [Ca/Fe], and 
[Ti/Fe]. In addition. Fig. [22] shows unweighted averages of those four element ratios where available. The agreement 
is good in all cases. Furthermore, the error bars seem to be reasonable estimates of the actual random and systematic 
error. 

The agreement between HRS and MRS [a/Fe] is very good (cr ~ 0.13 dex). Even though Fe lines outnumber a 
elements lines, the ratio [a/Fe] can be measured about as accurately as [Fe/H] because a and Fe respond similarly to 
errors in atmospheric parameters whereas Toff and [Fe/H] exhibit strong covariance. 

In addition to [Fe/H], B08b have published HRS measurements of [Ca/Fe]. Figure [23] shows the comparison between 
the stars we share in common (cr = 0.20 dex). The larger vertical scatter than horizontal scatter demonstrates that an 
MRS analysis is noisier than an HRS analysis when the number of measurable lines is small. Regardless, the degree of 
correlation is high, with a linear Pearson correlation coefficient of 0.5S, indicating that the medium-resolution spectra 
have significant power to constrain [Ca/Fe]. 
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Fig. 21. — Same as Fig. [T9l except for [Mg/Fe] (upper left), [Si/Fe] (upper right), [Ca/Fe] (lower left), and [Ti/Fe] (lower right). Only 
those measurements with estimated errors less than 0.45 dex are shown. 




Fig. 22. — Same as Fig. I19l except for an average of the a elements. 



Fig. 23. — Com parison between the HRS measurements of [Ca/Fe] 
of IBattaeha et a l. (2008b) and the synthesis-based medium resolu- 
tion measurements of [Ca/Fe] (this work) for the stars observed in 
both studies. 



TABLE 6 

Multi-Element Abundances in Sculptor 
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00'^59™21!6 


-33°42'58" 


19.662 ±0.090 


18.557 ±0.048 


5105 


2.02 


1.66 


-2.17 ±0. 


,19 




-0.12 + 0.75 


+0.39 + 0.28 




00'^59™21!7 


-33°41'02" 


18.296 ±0.037 


16.987 ±0.027 


4674 


1.27 
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-1.91 ±0. 
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+0.05 ±0.41 


+0.07 + 0.21 


+0.06 + 0.15 
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19.614 ±0.091 


18.383 ±0.048 


4805 


1.87 


1.70 


-1.95 ±0. 


,16 


+0.27 ±0.64 


-0.12 + 0.39 


+0.37 + 0.22 


-0.04 + 0.37 


001^ 59™ 26? 9 


-33°40'29" 


18.871 ±0.051 


17.559 ± 0.029 
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Note. — 



Table[6]is published in its entirety in the electronic edition of the Astrophysical Journal. A portion is shown here for guidance regarding form and content. 



